Exact Solution for Bulk-Edge Coupling in the Non-Abelian v = hl1 Quantum Hall Interferometer 
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It has been predicted that the phase sensitive part of the current through a non-abehan 1^ — 5/2 quantum Hall 
Fabry-Perot interferometer will depend on the number of localized charged e /4 quasiparticles (QPs) inside the 
interferometer cell. In the limit where all QPs are far from the edge, the leading contribution to the interference 
current is predicted to be absent if the number of enclosed QPs is odd and present otherwise, as a consequence 
of the non-abelian QP statistics. The situation is more complicated, however, if a localized QP is close enough 
to the boundary so that it can exchange a Majorana fermion with the edge via a tunneling process. Here, we 
derive an exact solution for the dependence of the interference current on the coupling strength for this tunneling 
process, and confirm a previous prediction that for sufficiently strong coupling, the localized QP is effectively 
incorporated in the edge and no longer affects the interference pattern. We confirm that the dimensionless 
coupling strength can be tuned by the source-drain voltage, and we find that not only does the magnitude of the 
even-odd effect change with the strength of bulk-edge coupling, but in addition, there is a universal shift in the 
interference phase as a function of coupling strength. Some implications for experiments are discussed at the 
end. 



I. INTRODUCTION 

Quantum mechanical systems with topological excitations 
are, in principle, ideally qualified for quantum information 
processing as the state of the topological sector is not acces- 
sible to local operators, and hence local perturbations cannot 
lead to decoherence. The best studied example of topologi- 
cally ordered states are fractional quantum Hall (QH) systems, 
where braiding of quasiparticles changes the ground state 
wave function. While in conventional QH states with odd de- 
nominator filling fraction only the phase of the wave function 
is changed by quasiparticle (QP) braids, the more recently dis- 
covered quantum Hall state at filling fraction 5/2 is expected 
to support non-abelian QPs, whose braiding corresponds to 
transformations in a degenerate ground state manifoldji^ The 
ground state degeneracy can be used to store information in 
the form of quantum bits, with one qubit for a pair of QPs, and 
quantum gate operations may be performed by the braiding of 
QPs. In principle, the read-out of quantum bits is possible by 
means of interference experiments.-^- 

One possible device for the readout of a topological quan- 
tum bit is the Fabry-Perot interferometer It consists of 
two narrow constrictions in a Hall bar which act as quan- 
tum point contacts and allow backscattering between counter- 
propagating edges modes. Interference between partial waves 
backscattered at the first and at the second quantum point con- 
tact is sensitive to the phase acquired during a trip around the 
interferometer cell. Proposals for using a Fabry-Perot interfer- 
ometer for the readout of a topological quantum bit rely on the 
fact that there is a relative change of the interference phase by 
TT depending on the state of the qubit enclosed in the cell,^^ 

The 5/2 quantum Hall state can be described as a p-wave su- 
perconductor of composite fermions.^ In this picture, QP exci- 
tations are vortices accompanied by an electric charge ±e/4. 
These vortices have a zero energy bound state at their core, 
which is described by a Majorana degree of freedom. The 



Majorana operators associated with two such vortices can be 
combined into a complex fermion and constitute a two level 
system suitable for quantum information processing. The 
state of this two level system can be changed by moving a 
third QP around one partner of the pair.- Depending on the 
occupancy of this two level system, the interference phase ob- 
tained by a partial wave that encircles the two level system 
is predicted to change by tt. More generally, the influence 
of bulk particles localized inside the interference cell on the 
Fabry-Perot interference phase can be used to provide evi- 
dence for the non-abelian character of 5/2-QP excitations. 

As an example, we may consider the phase change when 
voltage is applied to an ideal side gate, which is able to vary 
the area A enclosed by the interferometer path, without chang- 
ing the electron density inside. If the interference signal is 
caused by backscattering of e/4 quasiparticles at the con- 
strictions of the interferometer, and if there are no localized 
QPs inside the loop, then the phase of the interferometer sig- 
nal will change by 2it when the area is varied by an amount 
A A = 4<i>o / B, where B is the magnetic field strength and $o 
is the flux quantum for an electron. For an odd number of bulk 
QPs inside the interferometer cell, however, the leading sinu- 
soidal dependence of the interference current on this phase is 
expected to vanish, while it is restored for an even number of 
bulk QPsi^ii^ This dependence of the interference signal on 
the parity of bulk QPs in the interferometer cell constitutes 
the so-called even-odd effect. 

In order for the even-odd effect to be observable, it is nec- 
essary that the quantum state of the localized QPs remain in- 
dependent of time during the course of the current measure- 
ment. This can be a problem, in real systems, because of tun- 
nel coupling between the bulk QPs and the edge. In principle, 
one can imagine two types of tunnel couplings: tunneling of 
charged e /4 QPs into and out of the interferometer cellS, and 
coupling between bulk Majorana degrees of freedom and the 
Majorana mode along the edge. The former process should be 
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suppressed by the requirement of charge neutrality due to the 
Coulomb interaction, which is expected to be strong in small 
interference devices. The latter process is likely to be experi- 
mentally relevant, and it is this process which is the focus of 
the current paper. 

As tunnel matrix elements typically depend on distance in 
an exponential way, QPs localized near the sample edge can 
have a sizable tunnel coupling. In device geometries defined 
by etch trenches or top gates, the electron density typically 
changes from its maximum value to zero over a distance of 
many magnetic lengths, such that one can expect the filling 
factor to deviate from the ideal value of 5/2 in some region 
near the interfering edge states. However, any deviation from 
the ideal filling fraction 5/2 implies a finite density of QPs near 
the edge. Due to the spatial proximity to the interfering edge 
states, these localized QPs can have a significant tunnel cou- 
pling to the edge, and a realistic description of a non-abelian 
Fabry -Perot interferometer needs to take bulk-edge coupling 
into account. 

The case of one weakly coupled bulk QP was analyzed in 
perturbation theory by Overbosch and Wen^, and the case 
of two bulk QPs coupled to opposite edges was studied 
by the present authors.'^ While for weak coupling, an ana- 
lytic perturbative solution was possible, the strong coupling 
regime was analyzed numerically for a lattice model. It was 
found that at T = 0, the dimensionless parameter describ- 
ing the strength of the bulk-edge coupling may be written as 
hX^ /vn\e*V\, where A is the tunnel-coupling strength, is 
the velocity of the Majorana edge mode, V is the source-drain 
voltage, and e* = e/ 4 is the QP charge. For small values of 
this parameter, bulk QPs effectively decouple from the edge, 
while for a large coupling strength they are absorbed by the 
edge. In this manuscript, we present an exact solution for the 
influence of bulk-edge coupling on the magnitude and phase 
of interference in a non-abelian Fabry-Perot interferometer. 

For example, we analyze the case of a single bulk QP whose 
Majorana mode is coupled to one edge of the interferometer 
We find that the interference current can be exactly evaluated 
in the regime where the length b of the interferometer is short 
compared to hvn/e*V. The interference current is reduced 
relative to the case of no bulk impurity by the modulus of 
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and the interference phase is shifted by the phase of this ex- 
pression. Thus, in the presence of bulk-edge coupling the 
even-odd effect is modified in an important way: instead of 
the absence of the leading harmonic of the interference cur- 
rent, this harmonic grows with decreasing voltage at a rate 
which is enhanced relative to the behavior in the absence of 
a bulk Majorana mode (the intensity grows cx 1/V instead 
of oc 1/W in the high voltage regime). In addition, there 
is a universal phase shift of e*^/** in the low voltage regime 
relative to the high voltage regime. The result Eq. ([T]) is in 
agreement with previous perturbative and numerical results. 
The result Eq. ([T]i predicts interference contrast and phase for 
arbitrary coupling strength and should be relevant for the in- 
terpretation of = 5/2 interference experiments. 



The manuscript is organized as follows: in section II we 
introduce our model in a continuum formulation; in section 
III we discuss the implications of our exact solution for in- 
terference experiments; in section IV we describe the lattice 
formulation of the model; in section V we derive the exact 
solution of the lattice model; and in section VI we present an 
interpretation of this solution in terms of a resummed pertur- 
bation theory. In section VII we conclude with a discussion 
of our main results and some additional comments on their 
application to experiments. While sections I, II, III, and VII 
are intended for the more general reader, sections IV to VI are 
more mathematical and describe the derivation of our results 
in some detail. 



II. CONTINUUM MODEL 

The model we consider is a Hall bar parallel to the x-axis 
with two quantum point contactsi^i2ii2i]ii2ii2, see Fig. 1. In 
the absence of any coupling to bulk quasiparticles the upper 
(u) and lower (d) edges of the ly = 5/2 state are described by 
two charged boson fields <j>d{x) and two neutral Majo- 

rana fermion fields ipu{x), fpdix). The Lagrangian density for 
the boson field on each edge is that of a chiral Luttinger liquid 
with velocity Vc 



(2) 



Here, r = u,d denotes the upper and lower edge, and the 
plus sign goes with r = u, the minus sign with r — d. The 
charge density is given by pr = -^dxfr- For simplicity we 
set ?i = 1 when no confusion results. The Majorana fields 
encode the non-abelian properties of the 5/2-state, their La- 
grangian densities are 
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Within the p-wave superconductivity picture of the 5/2-state, 
QPs with a charge ±e/4 are vortices of the superconductor, 
which have Majorana bound states at their core. In our model, 
there are two localized bulk QPs which carry a zero mode 
Majorana each, described by a localized Majorana operator 
We denote the two bulk Majorana operators by r„, r^;, with 
the subscript indicating the edge to which the quasiparticle 
couples. In the absence of coupling, the Lagrangians of these 
bulk Majorana modes are 



(4) 



We assume the two bulk QPs to be spatially well separated 
from each other such that there is no coupling between the 
Majorana modes r„ and associated with them. The two- 
dimensional Hilbert space created by these Majorana modes is 
spanned by the two eigenvectors of the operator iTuTd- The 
coupling of Tu to the upper edge and to the lower edge, 
both at a; = a;o, is described by the Lagrangian density 



(5) 
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Bulk-edge coupling gives rise to tunneling times t\r = 
(7rA^/2z;„)^^. In order to judge the effect of bulk-edge cou- 
pling on the interference signal, the tunneling time has to be 
compared to the geometric time tf, = b/vn needed to move 
between the two constrictions separated by a distance b, and 
to the voltage time ty = fi/e*V, which can be interpreted 
as the extension in time of a QP wave packet. In the limit 
W ^ in which the interference signal is most clearly seen, 
the effective strength of bulk-edge coupling is given by the ra- 
tio tv/t\. As we shall see, if this ratio is much smaller than 
one, the bulk state is effectively decoupled from the edge, and 
if the ratio is much larger than one, the bulk state is effectively 
absorbed by the edge and does not influence the interference 
signal any more.-l^ 

A charge e /4 QP consists of a charge part and the neutral 
Majorana mode associated with it. The operator that tunnels 
a quasiparticle across a constriction is the product of a charge 
part and a neutral part which encodes the non-abelian statistics 
of quasiparticles. The tunneling part of the Hamiltonian is 



Ht^„ = f + f^, (6) 



where 



T : 



Ae'Vt 



VilClMl + mCRNRTTuTd] (7) 



transfers a quasiparticle from the lower to the upper edge 
through the left {L) and right [R) constrictions respectively, 
and its hermitian conjugate similarly transfers a quasipar- 
ticle from the upper to the lower edge. Alternatively, one can 
say that T transfers a quasihole from the upper to the lower 
edge, and that transfers a quasihole from lower to upper 
edge. The operators Ci.(fl) and AfL.(R) will be defined below. 
The Aharonov-Bohm phase is absorbed into the relative phase 
between the tunneling coefficients rji, and rjB- Here, V is the 
voltage difference between the two edges, and e* = e/4 is the 
quasiparticle charge. Correspondingly, the current operator is 
given by 



J=—{f-f^) 
i 



The operators 



(8) 



(9) 



are the charge part of the tunneling operator, operating on the 
charge mode. The factor 1/2 in the exponential reflects the 
fact that the QP charge e /4 is one half of the "natural" charge 
of a state with filling fraction v ^ 1/2. 

The neutral parts of the tunneling operators can be ex- 
pressed as spin fields of an Ising CFT'^ 



L(R) 



cj^{Tb/2)(Jd{Tb/2) 



(10) 



The (7 operators can be defined through their operation on the 
Majorana fermion fields asii 

il)r{y)cTr{xQ) = Tsgii(a;o - ?/)crr(xo)V'r(y) (11) 
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FIG. 1: Interferometer with quasi-particle tunneling at positions 
—b/2 and b/2. One localized Majorana mode couples to the up- 
per edge at spatial position x = xq, another one to the lower edge at 
the same position. For the derivation of the exact solution, the bulk 
Majoranas are positioned at the boundary of the interferometer cell 
with xo — b/2. 



with r = u,d, and the minus sign going with the upper edge. 
We will discuss an alternative expression for the neutral tun- 
neling operator in the next paragraph. The factor of iTuTd 
in the second term of Eq. ^ is included to account for the 
wrapping of a tunneling quasiparticle at position x = b/2 
around the two localized quasiparticles. This factor is respon- 
sible for the TT phase shift between the interference patterns 
corresponding to the two eigenvectors of iTu^d- 

The neutral part of the tunneling operator can be expressed 
in a more intuitive way by using the parity operator for the 
part of the system to the left of the tunneling site. We arrive 
at this formulation by representing the 5/2 state as a p-wave 
superconductor of composite fermions^. In this picture, the 
quasiparticle with charge e/4 is a vortex of the superconduc- 
tor The superconducting phase changes by 2tt when encir- 
cling the vortex so that the condensate wave function is single 
valued. As a Cooper pair has two fermions, the phase of the 
fermionic wave function changes only by tt when encircling 
the vortex, and there has to be a branch cut in the phase field 
seen by unpaired fermions to make their wave function single 
valued. For this reason, every vortex drags behind it a branch 
cut in the phase field of unpaired fermions. The phase jump 
of TT across the branch cut shows up as the minus sign in the 
commutation relation Eq. ( fTTl i: while a Majorana operator at 
spatial coordinate y > xq is not affected by the tunneling of 
a charge e/4 QP as described by a cr operator, for y < xq the 
Majorana operator acquires an extra minus sign. 

Alternatively, the minus sign a Majorana operator acquires 
when crossing the branch cut left behind by an e/4 QP can be 
described by an operator which shifts the phase of fermions by 
TT. This operator can be found by using an analogy with spatial 
translations. A spatial translation by a distance a is described 
by the exponential exp(ipa), where p is the momentum conju- 
gate to the spatial coordinate which is shifted by a. In order to 
describe a "phase translation", we need the exponential of the 
operator conjugate to the phase operator. In a superconductor, 
phase and number of cooper pairs are conjugate so that the 
Cooper pair number operator generates phase shifts. At zero 
temperature, the Cooper pair density is just half the electron 
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density p{x), and the operator 



ITT f ^ d^rp(r) 



(12) 



generates a relative phase shift of tt between operators ?/'(?/) 
with y < xo and y > xq, respectively. As N{—oo, xq) is the 
fermion number operator for the region x < xq, the opera- 
tor e*'^^'"""'^") is just the parity of the number of fermions 
to the left of xq, and we identify the operator P(— oo,a;o) 
defined in Eq. ( fT2b as the parity operator When evaluating 
P{—oo,xo)ipr{y) for y < xo, the change in particle num- 
ber due to the action of t/jriy) is included in the evaluation of 
P{—oo,xo) and changes its value by minus one, while for the 
opposite operator order ipr{y)P{-'Oo, xq) this is not the case. 
In this way, Eq. ( fTTT i is reproduced for y < xq when making 
the identification 



t(a;o)CTd(a;o) = P(-oo,a;o) 



(13) 



Clearly, for y > xq the order of P{—oo, xq) and t/jriy) does 
not influence the value of P(— oo, xq), and Eq. ( fTTI ) is repro- 
duced again. 

In the following, it will be useful to decompose the 
parity operator into a bulk part Pbuik measuring the par- 
ity of bulk Majoranas and an edge part Pcdgc- In or- 
der to keep our model simple, we assume that any local- 
ized QPs in the region x < —6/2 are far from the edge, 
so that the occupation number of their associated Majo- 
rana states do not change during the course of the exper- 
iment. Since the parity operator factorizes according to 
Pbuik (-00,6/2) = Pbuik(-oo,-6/2)Pbuik(-V2,V2), un- 
der the above assumption it is sufficient to include only 
-Pbuik(— fe/2, 6/2) in the tunneling operator for the right con- 
striction. In our model, there are only two bulk Majoranas Tu, 

inside the interferometer cell with —6/2 < x < 6/2. Their 
parity is determined by the operator iTu^d, which indeed ap- 
pears as a factor in the definition Eq. (|7]i of the tunneling op- 
erator 

As the bulk part of the parity operator is fully described 
by the factor iTu^d in the tunneling Hamiltonian, the neutral 
operators A/l,j? can be expressed in terms of the edge parity 
operator. In order to find an explicit expression for it, we ex- 
press the particle density on upper and lower edge together as 
Pcdgdx) = iil)u{x)'4)d{x) and find 
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(14) 



As the edge parity operator factorizes in the same way as the 
bulk parity operator, the equal time neutral correlation func- 
tion (AA(-6/2)AA(6/2)) = (Pedgo(-6/2,V2)) is given by 
the expectation value of the edge parity operator for the inter- 
ferometer cell 



Pcdgc(-6/2,6/2) 



(15) 



This expression will be useful for the treatment of the lattice 
model introduced in section IV. In the framework of this lattice 
model, the edge parity operator reduces to a product over local 
parity operators. 



III. SIGNATURES OF BULK-EDGE COUPLING IN THE 
INTERFERENCE CURRENT 

To lowest order in the tunnel couplings tir, the expec- 
tation value of the interference contribution to the backscat- 
tered current can be obtained using linear response theory. In 
this approach, the perturbation is the tunneling Hamiltonian 
Eq. ^ and O. Starting from the relation 



1 

ih 



dt{[J{0),Htunm 



(16) 



we find after some algebra that the interference contribution 
to the backscattered current is given by 



4e* 

^int = -^'R-eiT]lr]R 



[TrCl{T)CR{Q)) 



(17) 



where Tr is the time ordering operator, and 5 is the short time 
cutoff of the theory. Using the definition Eq. (|9]), the expecta- 
tion value of the charged correlator can be directly evaluated 
as 



(T,Ci(T)Cfl(0)) = 



^1/4 



62/^2] 1/8 



(18) 



In the absence of bulk-edge couphng, the neutral correlator 
can be obtained from the representation Eq. ( fTOl ) by using the 
the conformal dimension h„ ^ -^of the cr-field in the expres- 
sion for CFT correlation functionsi^. One finds 



{TrMl{r)AfRm 



51/4 



1/8 



(19) 



Alternatively, it can be calculated by using a bosonized ver- 
sion of Eq. (fT4] l. Using the lattice model described in the next 
section, we have been able to obtain an exact solution for the 
neutral correlation function in the presence of two impurities 
at a;o ~ |, one of them coupled to the upper edge and the 
second coupled to the lower edge. From Eq. dMI ) we find 



{TrNl{T)NR{Q)lTu{Q)Vdm 



2J4 



6^ 

„,2 



eAi(''-».-)^SXo(A2(6-^«„r)/^;2) 



(20) 



where A"o(a;) is the modified Bessel function of order zero. 
Using this expression in Eq. ( fTTI i. the interference current can 
be evaluated for arbitrary system size, bulk edge coupling, and 
all ratios of Vn/vc- For illustrative purposes, we evaluate now 
the interference current in the regime of small interferome- 
ter length 6, in which the interference contrast is highest and 
where the size 6 of the interferometer cell can be set to zero. 
In addition, we will concentrate on a situation with only one 
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impurity present in the bulk, say Tu- This situation can be de- 
scribed by sending oo in Eq. ( l20b such that the impu- 

rity degree of freedom Td is effectively absorbed by the edge. 
This situation has the benefit of being more easily interpreted 
than the two impurity case. We will see that the visibility of 
the interference signal grows from zero to one as the coupling 
strength is increased. We define oj = —e*V/h and find 



'int 



(21) 




-t'^ + i5t]' 
A?, 



with 



(22) 



1 



2x 



To evaluate the time integral in Eq. (|2T]) . we split the domain 
of integration into negative and positive times. In each sub- 
domain, the square roots of time in numerator and denom- 
inator cancel up to a phase factor so that the integral is a 
Fourier integral of Bessel functions, which can be found in 
the literatur&i^. For positive arguments, the Bessel function 
Kq{x) is real, while for negative arguments a; < 0, Ko{x) 
has a cut along the real axis and can be decomposed into real 
and imaginary part according to Ko{x — i6) = Ko{—x) + 
inloi-x). 

With respect to the case without buUc impurity, the interfer- 
ence signal is modified by the additional factor Jimp(A^/a;u„) 
in Eq. ( 1211 1. The modulus of Jimp reduces the amplitude of in- 
terference oscillations for small values of X^/ujVn, while the 
argument of Jimp gives rise to a phase shift. Both modulus 
and phase of Jimp{x) are displayed in Fig. 2. The expansion 
of Jimp (a;) for small arguments (weak tunneling or large volt- 
age) is 

Jimp(a;) = (l-^)^/^ + (1 + i) x^/^ + 0{xi) . (23) 

The expansion for large arguments (strong coupling expan- 
sion) is 



1 



Ax 



32x^ 



+ O 



(24) 



Interestingly, bulk edge coupling not only reduces the visi- 
bility of interference oscillations but also contributes a phase 
shift of 7r/4 between the weak and strong coupling limit. This 
universal phase shift as a function of voltage is a signature 
of bulk-edge coupling in a non-abelian interferometer. The 
phase e^*^/** characteristic of the weak coupling limit can be 
interpreted as the non-abelian part of the phase acquired by 
two charge — e/4 QHs encircling each other in the clockwise 
direction. This phase factor agrees with that obtained from the 
CFT correlation function of two cr-operators, which describe 




x-'h 
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FIG. 2: Modulus and argument of the factor Jimp(- __^ e*v ; 
ing the modification of the interference current due to bulk-edge cou- 
pling. 



the non-abelian part of charge ±e/4 particles. Alternatively, 
the non-abelian phase can be inferred from the fact that non- 
abelian QHs behave relative to each other as bosons if they are 
in the I fusion channel, i.e. that abelian and non-abelian part 
of the phase cancel each otheriiii^. Hence, the non-abelian 
phase has to compensate the abelian phase e''^/^ found from 
the operator product expansion of two e*'''^^)/^ QH operators. 



IV. LATTICE VERSION OF THE CONTINUUM MODEL 

In order to evaluate neutral correlation functions in the pres- 
ence of bulk-edge coupling beyond perturbation theory, we 
develop a lattice description of the continuum model intro- 
duced in section II. For the lattice model, the parity expecta- 
tion value can be evaluated numerically for arbitrary strength 
of bulk edge coupling'", and in section V we will derive the 
exact solution Eq. ( |20l i by using the inversion formula^ for 
so-called Hilbert type matrices. We shall first concentrate on 
the equal time correlation function 



(25) 



where P{—b/2, 6/2) is the edge parity operator for the up and 
down Majorana modes as defined in Eq. (fTSl l. As a first step 
towards defining a lattice version of this operator, we consider 
a one-dimensional model of complex lattice fermions defined 
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by the Hamiltonian 



(26) 



where the sum rans over lattice points j, and the lattice con- 
stant is denoted by a. The kinetic energy describes a disper- 
sion relation 



sufficient to consider the 7-Majorana mode only. We can now 
identify the right-moving branch of 7 with the upper edge and 
the left-moving branch with the lower edge. 

States with momenta between — tt and are occupied, 
hence the zero temperature correlation function is given by 



= 2i 



27r TT 



1 



I 



(31) 



6(fc) 



■ COS k 



'TT < k < TT 



(27) 



We study the model at half filling with kp = ^. Defining 
N = ^, the equal time parity expectation value for the com- 
plex lattice fermions is given by 



The correlation function vanishes if j — I is even, and is odd 
under exchange of j and /. The parity expectation value for a 
system of a right moving and a left moving Majorana mode is 
given by 



,-Eri!„/.c]c,^ ^ ^ TT (2ctc^._i)) . (28) (^b) = ( n (^^^0) = \/det((z7,7,'» -(32) 



It is expressed as the expectation value of a product of lattice 
operators. In order to make contact with a non-abelian 5/2- 
edge, we define Majorana operators 



7j 



. + g-.jV2 J 



2clci - 1 



(29a) 
(29b) 
(29c) 



This transformation corresponds to a boost to the right moving 
Fermi point. The left moving Fermi point now corresponds to 
the momentum — tt. Using the equality Eq. ( |29ct . the parity 
expectation value Eq. (1281 ) can be expressed as the expectation 
value of a product of 7 and 7 operators. As the Hamiltonian is 
quadratic, Wicks's theorem can be used to evaluate the expec- 
tation value of this product as the Pfaffian form of the matrix 
of correlation functions. The Hamiltonians for the 7- and 7- 
mode are 



n 



Aa 



4a 2n 
dk 

7_fc sin(fc) 7fc . (30b) 



2tt 



As there is no coupling term between the two Majo- 
rana modes, the correlation function matrix decomposes 
into a block with 7-correlators and another block with 7- 
correlators. Although we initially need both modes to write 
down an expression for the local parity operator, the deter- 
minant of the correlation function matrix factorizes and it is 



-<j<- 



The rh.s. of Eq. ( [32t is positive as required because the eigen- 
values of the matrix of correlation functions (ijj^fj') occur 
in pairs ±iA with real A. When evaluating the determinant 
for different system sizes numerically, on finds that it decays 
as N~^^^ in agreement with the analytical result obtained by 
using (T-correlators in the Ising CFT. 

We next want to calculate the influence of localized bulk 
modes on the parity expectation value. More specifically, we 
need to calculate the expectation value 



(33) 



-<j<' 



In order to evaluate this expectation value, we need to know 
the edge-edge, the impurity-impurity, and the impurity-edge 
correlation functions in the presence of a coupling between 
impurities and edge. The lattice version of the bulk-edge cou- 
pling Eq. ^ is 

^scatter = [ f{k)[i^Tulk + i d'^ -2k p +k\ ■ 

J 2tt ^ y/a y/a 

(34) 

Here, f{k) = f{—k) is unity for momenta fc ^ 1 and drops 
rapidly to zero for larger momenta, such that the dispersion 
relation can be linearized around the two Fermi points. As the 
Hamiltonian for bulk and edge Majorana states is quadratic, 
all correlation functions needed for the evaluation of the par- 
ity expectation value Eq. (|33] | can be evaluated exactly. After 
integrating out the edge Majorana modes, one obtains an ef- 
fective action for the bulk states 



T 



E 



' 2 



+ Af 



dq 

TT 



/((7)Go(ie„,g) r„(e„) +rd(-e„) 



'-Y + ^df ^/(g)Go(*6„, q ~ 2fc^^)V,(e„) 



(35) 



Here, 

Go(ie„,fc) = 



2 



denotes the edge Green function in the absence of scattering, 

(36) 
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e„ = (2n + l)kBT/h is a fermionic Matsubara frequency, 
and T denotes temperature. Although we use a finite tem- 
perature formalism here, we will focus on the zero temper- 
ature case in the end. As there is no coupling between the 
bulk Majorana states, the impurity-impurity correlator van- 
ishes {iVu^d) = 0. The correlation functions between im- 
purity and edge operators are given by 



(ir„(-e;)7fc(e,)) = 



(ird(-e;)7fc(ei)) 



x(r„(-eOr„(eO>5.„ 

Go{iei,k-2kF) f{k) ^ 

x(rrf(-Q)rd(Q)><s„ 



(37) 



Due to chirality, the edge-edge correlation function {ijiji') 
depends on the bulk-edge coupling strengths A^, A^ only if 
the coupling to the impurity occurs between the two lattice 
sites I and I'. As only edge-edge correlation functions with 
both I and I' inside the interferometer cell are needed for the 
evaluation of Eq. ( |33] ), these edge-edge correlators become in- 
dependent of the bulk-edge coupling strengths for xq = 6/2, 
i.e. for impurities coupling to the boundary of the interferom- 
eter cell. In order to simplify the task of calculating the full 
neutral correlation function including bulk-edge coupling, we 
will adopt a;o = 5/2 in the following. 

In order to extract the universal long distance behavior 
of correlation functions, we linearize the dispersion relation 
around the two Fermi points and remove the momentum cut- 
off when possible. In this way, the Fourier transform fj of a 
function / sin fc] becomes 



dk 



37r/2 



27r 



— sm k 
a 



(- 



1 



dk 
2^ 



/ 



v„ 



(38) 



The first term with a rapidly oscillating position dependence is 
due to integrating over momenta — 37r/2 < k — 7t/2, whereas 
the second term with a smooth position dependence is due 
to integration over momenta — 7r/2 < k < tt/2. With the 
help of this formula and regularizing momentum integrals as 
^ Jg°° dfc e*'^'-^^*''-' = ^ x+iTj ^ expression Eq. (ISTT l is easily 
reproduced in the limit ry ^ 0. For the bulk-edge correlation 
functions, one finds 



(39) 



(zr,7.) = (-1)^- ^^e^^^^"/""i?i(2A^,a/.,2j.(40) 



Here, the exponential integral is defined as 



noo p — tx 

Ei{x) = / dt 



t 



(41) 



It has the asymptotic expansions 

Ei{x) -7-lna; for a; ^ , (42) 

El (x) for a; — > oo 

X 



with 7 = 0.57721... denoting Euler's constant. Although it is 
not needed for the calculations presented in this manuscript, 
we would like to mention the result for the edge-edge corre- 
lation function in the presence of bulk-edge coupling, which 
was used to obtain the numerical result for the reduction factor 
in Ref. [13 It is given by 



1 1 



(43) 



AXla 



V. 



X mmi-i) - e(-j)e(0] 

-(-ly-' ^ e2^S('-^>/-" Ei[2\l{l^j)a/vl] 

X M-mi) - e(j)e(-/)] . 



EXACT SOLUTION FOR PARITY CORRELATION 
FUNCTION WITH IMPURITIES 



As explained in the last section, we consider a geometry 
where the bulk impurity couples to the edge at the boundary 
of the interferometer cell, i.e. a-o = h/2. This geometry some- 
what simplifies calculations because now the edge-edge cor- 
relation function does not have an impurity contribution. For 

ease of notation, we assume = = A^ in the follow- 
ing, the generalization to two different couplings is straight- 
forward. 

Again we will be calculating a correlation function by using 
Wick's theorem to rewrite that correlation function as a deter- 
minant analogous to Eq. ( |32] |. However, here we will calcu- 
late the more complicated correlation function Eq. ( [33] ). We 
denote the matrix of correlation functions, whose determinant 
needs to be calculated, by C. All diagonal elements of C van- 
ish. We adopt a bra-ket notation in the following and denote a 
position along the edge by \ and the two bulk impurities by 
\u), \d). The edge-edge correlation function is then given by 
(j|C|j'), the bulk-edge correlation is 



O'lci^.) 



2A 



Eii2X^j) 



2A 



{j\C\d) = i-iy-e^^ 'Ei{2X^j) 



(44) 



The impurity-impurity correlation function is (m|C|c?) = 
((i|C|u) = 0. Our result is again a square root of a determi- 
nant, and the determinant is the product of all eigenvalues. Let 
us assume we know the eigenvalues of the edge-edge part. For 
small A, the bulk-edge correlation functions (j|C|u), {j\C\d) 
are of order A In^ A, and the leading contribution to the deter- 
minant is obtained by multiplying the determinant of the edge- 
edge correlators with the perturbatively calculated eigenval- 
ues of the impurity-impurity part of the matrix. Although the 
eigenvalues are calculated perturbatively in A, which is pro- 
portional to the lattice constant and goes to zero in the contin- 
uum limit, the final result is valid even in the strong coupling 
regime with large AA^. The square root of the product of these 
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two eigenvalues is the reduction factor R{b), i.e. the ratio 



R 



(45) 



of the neutral expectation value in the presence of two impu- 
rities to the expectation value without impurities. 

Without bulk-edge coupling, C has two zero eigenvalues. 
To determine the shift of these zero eigenvalues due to the 
coupling between impurities and edges, we use second order 
perturbation theory to calculate the effective matrix elements 
(w|C|(i)off due to "virtual transitions" of a bulk Majorana to 
the edge and back. Up to a sign, the reduction factor is then 
equal to this effective matrix element. 



R 



{u\Q\d), 



(46) 



To calculate the effective matrix element, we change to a new 
basis 



\a) = \u) + \d) , \b) = \u) - \d) 



(47) 



In the new basis, \a) only couples to even lattice sites, while 
\b) only couples to odd ones. To exploit this, it is useful to 
decompose the lattice into even and odd sites according to 

J = 2n + i(l + a) , / = 2n' + ^(1 + a') . (48) 

We assume that the number of lattice sites N is even such that 
this decomposition works. Then, n runs from to N/2, and 
cr = ±1 determines whether the lattice site is even or odd. In 
analogy to Eq. ( |48l ), we use an \na) basis in the following, 
where |ncr) = \j) with j given by Eq. ( l48T l. Denoting the 
eigenvectors of the bulk-bulk part of C by |e;) and the 
corresponding eigenvalues by A/, the effective matrix element 
between \a) and \b) state is given by 



N 



{b\C\aU = ^ 



1=1 



(6|C|e;)(e;|C|a) 
Az 



(49) 



= J2(^\C\j){j\{C"r'\f){f\C\a) 

= Y.mn+){n+\{C)-'\n'-){n' -\Q.\a) 

nn' 

with (see Eq. ( l40l i 

(6|C|n+) = (n - |C|a) ^ —e^^'^^E^ UX^n) . (50) 

The physical interpretation of the effective matrix element 
Eq. ( |49] ) is that the \a) bulk state makes a virtual transition 
to the edge and then back to the \b) state. For this reason, the 
double sum in the second line of Eq. ( |49l ) runs over edge states 
only. In order to calculate matrix elements of (C°)^^, we use 
the fact that in the n, a basis C° has the form 



C0 = 



D 

-D^ 



(51) 



with 



D 2 1 



see Eq. dlTl i. Then, 



|C-i|n'-) = (D-n 



(52) 



(53) 



Matrices of the form of D are known as Hilbert-type, and us- 
ing the inversion formula derived by Trench Scheinobi^ we 
find that the inverse of D is given by 



(D-i) 



TT 

4 n — m 



'—jTl 



1 



2 q^rn 

1 

2{n - 



1 



1 



2(g - m) 



(54) 



1 1 r(i + m)r(i + f 

Trn - m - i T{m)V (l + f 

rfn-i)rf| + ^-n) 



r(n)r 1 



N 



n) 



(55) 



As we will finally take the limit N ^ oo with A^A^ fixed, we 
can use Sterling's formula to simplify 



1 



1 



TT n — m 




1 + f -m 



(56) 



As the formula Eq. ( |49] l uses only the symmetric part of D ^, 
we symmetrize and obtain 



N + 2 



1 



even,mn 



4tt 



mnJ(l + f -n)(f + l-m 



(57) 



Taking everything together, the reduction factor is 



^^A^/4Ay 

4tT \ TT J 



Tl-l 



N 



1-71 

(58)" 



Taking the continuum limit, one finds 



R 



X-'N— 
2 r 



/ dxe^'^Ei {2x) 
Jo 



^xiX^N ~x) 



X^N- 



TT 



Ko (VN) 



(59) 



For the evaluation of the integral in the last equation, we used 
the integral representation Eq. (HTI) for Ei (x) and evaluated 
the .T-integral in terms of a modified Bessel function Iq (x), for 
details see Ref. [l^ The remaining integral is again tabulated 
in Ref. [Ill One sees that the reduction factor is the square of 
reduction factors due to the two individual impurities. Using 
the asymptotic behavior of the zeroth order modified Bessel 
function Kq{z) w — Inz for z and Ko{z) w V?"^"^ 
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for z we find the asymptotic behavior of the reduction 

factor 



i?(A^iV<l) = -X^N{lnX^Ny 
R{X^N ^ oo) 



TT 
1 



(60) 



Because of the chirality of upper and lower edge, the time 
dependence of the exact solution can be obtained by replacing 
b ^ 6 — ir in the factor describing the upper edge and b 
b + ir in the factor describing the lower edge in Eq. ( |59] l such 
that the time dependent reduction factor is given by 



Rib,T) 



(61) 



- \/b'^ + v'^T^ 

'e>■>+'--^y<K,{\l{b + ^v^T)/vl 
'e'l^'--^y<K,{\l{b~^v^r,vl) 



This extension of the static solution to finite time differences 
can be justified by considering the case of one bulk impurity 
coupled to one of the edges, say the upper one, first. Then, we 
define a neutral correlation function 

(62) 

which in principle could depend on the four arguments xi, X2, 
Ti, T2 and the parameter xq (position of bulk-edge coupling) 
separately. As the impurity is static, there is translational in- 
variance in time and Mu can only depend on the difference 
T2 — Ti. As we consider a situation where the impurity is cou- 
pled to a point inside the cell, we have X2 < xq. Then, due 
to chirality, the field (7u{x2,T2) is not influenced by bulk-edge 
coupling and does not depend on X2 and T2 separately but only 
on the combination X2 + iT2, and it satisfies the differential 
equation {dx2 + idr2)<^u{x2, T2) = 0. If we restrict ourselves 
to time differences r2 — ti 7^ different from zero, the corre- 
lation function Mu satisfies the same differential equation and 
can for this reason only depend on the variables X2+i{T2—Ti) 
and xi. However, from Eq. ( l59l ) we see that the static corre- 
lation function only depends on the difference xi — X2, so we 
can conclude that Afu{x2, T2; a;i, n; xo) is a function of the 
single variable xi — X2 — «(t2 — ti), and that the correct an- 
alytic continuation of Eq. ( |59l ) is indeed given by replacing 
b ^ b — IT with b ^ xi — X2 and r = r2 — ti. The analytic 
continuation for the lower edge can be derived by a similar 
argument. 



INTERPRETATION IN TERMS OF RESUMMED 
PERTURBATION THEORY 



VI. 



In this section, we show that the exact solution Eq. 
can be reproduced by resumming the perturbative expansion 
of the neutral correlation function in powers of the bulk-edge 
coupling constant. The terms contributing to this resumma- 
tion are those which turn the zeroth order bulk-bulk correla- 
tion function (r^rr(T)rr(0))o ^ 1 into the full correlator. 



The (T„cr„^„-correlation function appearing in the lowest or- 
der expression Eq. (|64] | is not modified in the perturbative ex- 
pansion due to the special choice xq = 6/2, which implies 
that the chiral ctuO-mt/'u -correlation function is only evaluated 
for spatial arguments x < xq. As one can see for example 
from the edge-Majorana correlation function Eq. (l44l i. bulk- 
edge coupling is only important if one spatial argument is to 
the left and another one to the right of xq. 

Since upper and lower edge decouple in perturbation theory 
(modulo fusion channels), we consider only one impurity cou- 
pled to one edge, say the upper edge. We start by calculating 
the neutral equal time correlation function 

{^bl2)c7u{bl2)ViTu) (63) 
in perturbation theory. The lowest order contribution isiS 

/oo b b b 

dT{Tra.4--,0+)au{-,0~)M^.r))o 



x(T,r„(T)r„(0))o 



(64) 



This expression is logarithmically divergent and needs to be 
regularized by a cutoff on the time integral, which in Ref. [l^ 
was inserted by hand. However, when resumming the in- 
finite set of diagrams which turns the zero order correlator 
into the full correlator, the integral becomes finite and we are 
able to reproduce the exact solution for the reduction factor in 
Eq. i5% . In order to verify this proposition, we calculate the 
expression Eq. ( l64l ) with (T'^r„(T)r„(0))o replaced by the 
full correlation function (TT-r„(r)r„(0)). From the action 
Eq. (I35] ) we find in frequency space 



(65) 



-signej 



After calculating the Fourier transform to Matsubara time we 
find 



(T,r„(T)r„(o)) 



1 



-Im 



(66) 

The calculation can be expressed in a more compact fashion 
by defining 



1 



g{-z) = -e-2.^"/-"i?l(-2zA2/^;„) 



TT 

In addition, we make use of the CFT correlation function 

1 



(67) 



V2, 0+)a„(6/2, 0-)^'„(6/2, r))o 



/27r 



1 



y/{b + iT)iT 



(68) 



Note that the additional factor 1 / in Eq. (I68b as compared 
to Refill is due to the difference in the Majorana Lagrangian 
Eq. dS]) as compared to\l^ Now we can express the equal time 
correlation function as 



lu = \fi\u{-b)i 



1 



df_ g{-z)-g{z) 
i ^yz{z + b) 



(69) 
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We close the integration contour in the right half plane for the 
integral over g{z) and in the left half plane for g(—z). We are 
allowed to do so because g{z) ~ 1/z for large z in the right 
half plane, which together with the asymptotic 1/z behavior 
of the (criT^') correlator makes sure that the infinite semicircle 
does not contribute to the integral. As the (crcr^') correlator 
has a cut only along the negative real axis between z = —b 
and z = 0, the integral over g{z) vanishes. The integral over 
g{—z) can be converted into a contour encircling this cut and 
gives 



ax : 

y'{b-x)x 



(70) 



We note that the expression for Td differs from Eq. (iTOl i by 
a factor of i due to the opposite chirality in the correlator 
Eq. ( l68l l. For this reason, the reduction factor obtained from 
the product XuXd is real in agreement with Eq. ( |59l ). The ex- 
pression Eq. dTOl i agrees up to a phase factor with the corre- 
lation function obtained from the square root of the reduction 
factor Eq. ( l59l ). As an additional benefit, using the finite tem- 
perature expressions for the (crcr^') and (r„rti) correlators 
opens a route towards a generalization to finite temperature. 



VII. DISCUSSION 

In this paper, we have analyzed the influence of a tunnel 
coupling between bulk and edge Majorana states on the vis- 
ibility and phase of interference oscillations in a non-abelian 
ly = 5/2 quantum Hall interferometer. Such a tunnel coupling 
is important because it blurs the distinction between bulk and 
edge degrees of freedom and thus complicates the observation 
of the even-odd effect as a signature for non-abelian statistics. 
In our discussion, we have focused on the behavior at temper- 
ature T = 0, for an interferometer encircling one or two lo- 
calized quasiparticles (QPs), as a function of the source-drain 
voltage and the strength of coupling between the Majorana 
modes of the impurities and the neutral modes of the edge. 

The present paper is an extension of results presented in a 
previous letter by the authors jiS, In the present paper, we have 
found an exact analytic formula for the equal-time parity cor- 
relation function for the two ends of the interferometer, when 
the localized quasiparticles are both located close to one end, 
and we have verified that the correlation function saturates, 
in the strong coupling limit, at the same value as one would 
find in the absence of localized quasiparticles. This correla- 
tion function was only obtained numerically in our previous 
work. Analyzing the analytic properties of the function in 
the space-time plane, we now obtain an analytic form for the 
correlation function at two different times, and from that we 
can predict the dependence of the interference amplitude on 
the applied voltage and the bulk-edge coupling strengths. We 
have also been able to examine the phase shift in the interfer- 
ence pattern introduced by the presence of a finite coupling 
between the bulk QPs and the edge. In addition, the current 



paper presents various details of the analysis that had to be 
omitted from Ref. [T^ due to lack of space. 

In our work, we have particularly examined the case of a 
short interferometer, or relatively low voltage, where the inter- 
ference visibility is largest. Specifically, we assume ty ^ h, 
where, = b/vn is the time needed for a neutral excitation 
to move along one edge, from one constriction the other, and 
tv — h/e*V is the extension in time of a QP wave packet 
transferred from one edge to the other by backscattering at 
one of the constrictions. We shall summarize here the quali- 
tative results of our studies, and then say a few words about 
their implications for experiments. 

We start by discussing the case of a single bulk Majorana 
mode inside the interferometer. In the absence of bulk-edge 
coupling, the leading harmonic of the interference signal van- 
ishes. In agreement with previous analyses of this problero^jiS, 
we find that for weak coupling, at T = 0, interference can be 
observed but is reduced by a factor proportional to \JTvJt\, 
where t\ = (7rA^/2w,i)^^ is the characteristic tunneling time 
associated with the exchange of a Majorana particle between 
the localized QP and the edge. For large values of the ef- 
fective coupling constant ty /t\, the bulk Majorana mode is 
effectively absorbed by the edge and the interference signal is 
fully restored, such that the strong coupling case corresponds 
to an interferometer with no bulk degree of freedom. In ad- 
dition, on the way from weak to strong coupling, the phase 
of the interference signal is shifted by it/ A. Although bulk- 
edge coupling enforces a modification of the way one looks 
at the even-odd effect, it actually enriches this effect with a 
new direction in parameter space, as the dimensionless cou- 
pling strength ty / 1 a depends on source drain voltage. The 
signature for an odd number of impurities, with just one of 
them coupled to the edge, is not the complete absence of the 
leading harmonic, but rather a reduced amplitude, which de- 
pends on the applied voltage. The interference intensity grows 
with decreasing voltage at a rate which is enhanced relative 
to the behavior in the absence of a bulk Majorana mode (the 
intensity grows cx 1/V instead of cx l/\/V ) until the reduc- 
tion factor saturates at unity. In addition, when the reduction 
factor is small compared to unity, the interference factor will 
have universal phase shift of 7r/4 relative to the pattern in the 
strongly coupled regime. 

Our results for a single bulk QP at T = can be readily 
generalized to the case of finite temperatures. Qualitatively, 
the temperature will make relatively little difference as long as 
kT is small compared to e*V . However in the opposite limit, 
temperature will be important, and its effect may be roughly 
estimated by replacing the the voltage time ty, in the formulas 
above, by the thermal time 1^ = H/kT. A more quantitative 
analysis can be given, but it will not be discussed here. 

In the case of two bulk Majorana modes coupled to opposite 
edges, with comparable coupling strengths, the average value 
of the interference signal shows the same qualitative behavior 
as in the case of a single impurity. However, the average value 
of the interference current now has a more subtle interpreta- 
tion than in the case of a single impurity: due to the presence 
of the factor iTu^d in the tunneling operator, the interference 
signal is sensitive to the state of the two bulk Majorana modes 
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even in the absence of bulk-edge coupling. Although in the 
absence of bulk-edge coupling the expectation value of iT^Td 
does not change in time once it is prepared in an eigenstate, 
the quantum statistical average {iTurd) is zero, as the average 
is taken over both possible states of the system. Experimen- 
tally, then, the quantum statistical average corresponds to a 
situation where the interference signal is averaged over differ- 
ent initializations of the bulk states. 

In the presence of weak but non-zero coupling between the 
bulk and edge modes, two things happen: i) the quantum me- 
chanical average (iT-uFd) oc (ty/tx) In^ [tx/ty] is now finite; 
and ii) for a single experiment starting with the bulk states ini- 
tialized in an eigenstate of iTu^d, the interference phase will 
fluctuate on the time scale mm{t\u, t\d), so that the time aver- 
age of the interference current over these quantum mechanical 
fluctuations is equal to the quantum statistical average. 

Now we must distinguish between two experimental situ- 
ations. If the experimental measurement is averaged over a 
time tav that is long compared to one or both of the switching 
times tx, then the experiment will measure the statistical av- 
erage interference signal, which will be very close to zero if 
t\ is much longer than ty and <t- On the other hand, if the 
bulk modes are so weakly coupled to the edge that for both 
bulk QPs, t\ > tav, then the experiment will not measure a 
statistical average, but will see only one or the other of two 
possible fermion states. In this case one will see a full inter- 
ference signal, with the same amplitude as if the impurities 
were not there at all, but with a phase that depends on the 
starting configuration of the system. 

More generally, we may distinguish three ranges of cou- 
pling strengths for a QP localized in the bulk. If coupling 
to the edge is so weak that tx > ti^-^, we may say that the 
bulk state is "decoupled" from the edge. If the coupling is 
in the range where t\ is small compared to ty and tx, then 
QP is "strongly coupled" to the edge. If the coupling strength 
is in the range where tx is small compared to tav, but large 
compared to either ty or tr, then we may say that the QP is 
"weakly coupled" to edge. 

In an experiment where the interfering edge encloses two 
or more localized QPs, we may ignore any strongly coupled 
QPs, as they will be effectively incorporated into the edge. If 
there are no "weakly coupled" QPs inside the loop, then the 
presence or absence of the interference signal corresponding 
to charge e/4 will be determined by the number of decoupled 
QPs inside the loop. The interference pattern will be present 
if this number is even, and absent otherwise. If there are any 
weakly coupled QPs inside the loop, however, with relaxation 
time tx small compared to the experimental averaging time, 
then the interference pattern will be absent, regardless of the 
number of decoupled QPs that may also be present. 

In a recent experiment by Willett et ali^, resistance oscilla- 
tions in a Fabry-Perot device were studied experimentally. For 
magnetic fields near a bulk filling fraction v = 5/2, oscilla- 
tions in the longitudinal resistance were observed as a function 
of side gate voltage. Depending on the range over which the 
side gate voltage was varied, consecutive doubling and halv- 
ing of the voltage period of resistance oscillations was ob- 
served. The side gate voltage was interpreted as changing the 



area of the interferometer cell. If the bulk filling fraction devi- 
ates slightly from the exact value of 5/2, a change of area will 
once in a while change the number of QPs inside the interfer- 
ometer cell by one. If the number of QPs changes from even 
to odd, the fundamental harmonic is suppressed and the period 
of the interference signal is halved, while for a change from 
an odd to an even number of localized QPs inside the cell, the 
period is doubled. A reduced voltage period in the presence 
of an odd number of localized QPs can arise from interference 
of abelian charge ±e/2 QPs or from non-abelian charge ±e/4 
QPs encircling the interferometer cell twice. This interpreta- 
tion of the experiments^ is discussed in more detail in Ref . lU. 

A difficulty with this interpretation is that it seems to re- 
quire that all QPs inside the edge are either completely de- 
coupled from the edge, or strongly coupled so that they are 
incorporated into the edge. If there any QPs in the intermedi- 
ate weakly coupled range, where the Majorana state changes 
back and forth many times during the averaging time of the 
experiment, then the interference corresponding to charge e/4 
would be completely absent, or at least greatly reduced in size. 
It is not clear to us, why there should be no weakly coupled 
QPs inside the interferometer in the interference experiments 
of Ref. H^, nor is it clear why one should be rid of their effects 
if such QPs are present. We do note, however, that one pos- 
sible ingredient that is missing from our analysis is a tunnel 
coupling that allows for an exchange of Majorana fermions 
between bulk quasi-particles. In the presence of such a cou- 
pling, the degeneracy of the localized bulk Majorana states is 
removed. The resulting spectrum is then composed of sev- 
eral states. Each of these states corresponds to an interfer- 
ence pattern, whose amplitude and phase are determined by 
the expectation value and fluctuations of the parity operator 
in that state. The coupling of the bulk states to the propagat- 
ing Majorana mode of the edge introduces a width to these 
states, equilibrates them to the temperature of the edge, and 
thermally averages the interference pattern. Then, if the split- 
ting between the bulk states is large compared to their width 
and to the temperature, a well defined interference pattern will 
be observed. A recent estimate?^ for the tunnel splitting be- 
tween quasi-particles is rather sizeable, about lOOmJ-iT for a 
separation of 0.1 micron. This effect and the overall problem 
clearly require further investigation. 

Altogether, the experimental results obtained so far are not 
yet well understood by theory, and it is not yet clear that the 
observed effects originate from the unique properties of non- 
abelian quasi-particles at = 5/2. It is possible that a clue to 
that question may be obtained from the transition region be- 
tween two different periods, when a change in the number of 
localized QPs happens due to a change in area of the interfer- 
ometer cell, and when the distance between edge and bulk QP 
decreases as a function of side gate voltage. Then, the bulk- 
edge coupling should change from small to large values, and 
one can expect that the theory developed in this manuscript 
is applicable. It would be interesting to analyze the data of 
Ref. from the point of view of bulk-edge coupling, and to 
test the theoretical prediction that in the transition region be- 
tween different gate voltage periods the interference current 
has a modified power law dependence on source-drain voltage 
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and that there is a 7r/4 phase shift as a function of source-drain 
voltage. 
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